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Abstract 

This paper reviews the importance of numerical drag 
prediction in an aircraft design environment. A 
chronicle of collaborations between the authors and 
colleagues is discussed. This retrospective provides 
a road-map which illustrates some of the actions 
taken in the past seven years in pursuit of accu- 
rate drag prediction. The advances made possible 
through these collaborations have changed the man- 
ner in which business is conducted during the design 
of all-new aircraft. 

The subject of this study is the DLR-F4 
wing/body transonic model. Specifically, the work 
conducted herein was in support of the I st CFD Drag 
Prediction Workshop, which was held in conjunction 
with the 19 t/l Applied Aerodynamics Conference in 
Anaheim, CA during June, 2001. 

Comprehensive sets of OVERFLOW simulations 
were independently performed by several users on 
a variety of computational platforms. CFL3D was 
used on a limited basis for additional comparison on 
the same overset mesh. Drag polars based on this 
database were constructed with a CFD-to-Test cor- 
rection applied and compared with test data from 
three facilities. These comparisons show that the 
predicted drag polars fall inside the scatter band of 
the test data, at least for pre-buffet conditions. This 
places the corrected drag levels within 1% of the 
averaged experimental values. At the design point, 
the OVERFLOW and CFL3D drag predictions are 
within 1-2% of each other. In addition, drag-rise 
characteristics and a boundary of drag-divergence 
Mach number are presented. 

*AIAA Associate Fellow, DPW Organizing Committee, 
Boeing Technical Fellow 

*AIAA Associate Fellow, Research Scientist 

* AIAA Associate Fellow, Research Scientist 
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Nomenclature 
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Wing Aspect Ftatio = ~ — 
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Acoustic Speed 
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Wing Span 

Cd 

Drag Coefficient = Dr °-2 
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C l 

Lift Coefficient = — ' t 

Cref 

Wing Reference Chord 

count 

Drag Coefficient Unit = 0.0001 

D 

Drag 

DPW 

Drag Prediction Workshop 

e 

Oswald’s Efficiency Factor 

L 

Lift 

M 

Mach Number 

RANS 

Reynolds- Averaged Navier-Stokes 

Ren 

Reynolds number = — V J° ° re S 
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Wing Reference Area 

SFC 

Specific Fuel Consumption 

W 

Weight 


Wall Distance — ^ ■ u y 

fl w 

Q 

Dynamic Pressure = |pV 2 

a 

Angle of Attack 

Ac/4 

Wing Quarter-Chord Sweep 

oo 

Signifies Freestream Conditions 


1 Introduction 

In reflection, as the 100-year anniversary of flight 
draws near, it is truly amazing just how far the in- 
dustry has progressed. This rapid advance in the 
science, technology, and business of flight was made 
possible through a blend of competition and coop- 
eration between industry rivals, government agen- 
cies, and academic institutions around the world. It 
is in this spirit that the Drag Prediction Workshop 
(DPW) was conducted [l]-[2]. Participants from six 
nations came together for a common goal - to assess 
the state-of-the-art of drag prediction using Compu- 
tational Fluid Dynamics (CFD) methods based on 
the Reynolds- Averaged Navier-Stokes (RANS) equa- 
tions. 
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Why is drag prediction so important to the indus- 
try of flight? In short, drag represents the irrecov- 
erable aerodynamic losses associated with a flight- 
based mission. For example, consider the generic 
task of delivering a payload between distant city 
pairs. The Breguet-Range equation, which aptly ap- 
plies to long-range missions of jet aircraft, is: 


Range — 


ML a 
D SFC n 


( W 0 + Wf \ 

V Wo ) ' 


(1) 


Here, M is the cruise Mach number, L & D are the 
aerodynamic forces of lift and drag, respectively, a 
is the acoustic speed, SFC is the specific fuel con- 
sumption of the engines, Wo is the aircraft landing 
weight, and Wj is the weight of fuel burned during 
the flight. The Breguet-Range equation illustrates 
the importance of drag prediction as a function of 
lift and Mach number in the context of aerodynamic 
design; it also provides a glimpse into the interplay 
between the various disciplines. 

Referring to Eqn (1), one might assume that the 
aerodynamic efficiency of an aircraft is represented 
by the propulsion efficiency is embedded in 

SFC , and that the structural efficiency directly im- 
pacts Wq . Interestingly, historical trends of in- 
service transport aircraft indicate that very little im- 
provement in the metric has been accomplished 
in the past 50 years. Yet it would be somewhat 
naive to state that no aerodynamic advances have 
been made during this period. In actuality, improve- 
ments in aerodynamics have better served aircraft 
designs by trading them for improvements in other 
disciplines. For example, the ability to increase the 
thickness-to-chord ratio of a wing while maintain- 
ing not only reduces the structural weight of 
the wing, it also provides additional fuel volume. 
In terms of Eqn (1), an aerodynamic improvement 
of this nature would manifest itself as a decrease in 
W 0 and an increase in Wj with the net result being 
an increase in range. Reducing the aircraft’s empty 
weight has the added benefit of reducing the cost of 
the vehicle. Obviously, this aerodynamic improve- 
ment would not be apparent in the trend charts of 
ML 
D ■ 

Assume that an airline would like to provide a ser- 
vice between two cities with an aircraft that, when 
fully loaded with payload and fuel, is 1% short on 
range. Since the aircraft is fuel-volume limited, the 
only recourse is to reduce the payload weight. In 
relative terms, a typical ratio of weights might have 
^7 3 W 0 and Wp a yi oa d — §Wo. In this scenario, 

Eqn (1) shows that the operator would have to re- 
duce the payload (read revenue) by 7.6% to recover 
the 1% shortfall on range. Since most airlines op- 
erate on very small margins, this service most likely 
will no longer be a profit-generating venture. This 


example illustrates that in the current business of 
flight, a 1% delta in aircraft performance is a sig- 
nificant change. While improving an aircraft’s per- 
formance by 1% may not be a trivial task given the 
usual constraints, losing 5% is easily done if attention 
is not paid to detail (e.g., juncture flows, external 
doublers, gaps, etc.). 

Now consider a more typical case where the air- 
craft does not suffer from a shortfall on range. In 
round numbers, the Direct Operating Costs (DOC) 
of a transport aircraft can be itemized as: 50% for 
the cost of ownership, 20% for fuel burn, 20% for 
crew salaries and maintanence, and 10% for miscel- 
laneous other items. From an airline’s perspective, if 
the DOC of its fleet of aircraft could be reduced by 
5% with a new design (while providing the same set 
of services to its customers), the airline would most 
likely retire its entire fleet and replace it with the 
new aircraft [3]. 

So how can aerodynamics be leveraged to im- 
prove the economics associated with a flight-based 
mission? A simplified high-lift-system design that 
retains and Cl max reduces manufacturing and 
maintanence costs as well as part count. Increas- 
ing the cruise Mach number without reducing pr 
reduces the time-dependent costs such as crew and 
maintanence. And the classic, increasing with- 
out penalizing the other disciplines reduces fuel burn. 
These are just a few examples of how aerodynamic 
advances can have an impact on DOC., .and all of 
these require accurate drag predictions. 

To push aerodynamic technologies forward, it is 
becoming more important that accurate drag predic- 
tion become a consistent product of the CFD com- 
munity. Once this prerequisite is accomplished, the 
full benefits of automated aerodynamic shape opti- 
mization may begin to be realized. 

With the various on-going design programs, these 
are exciting times for the aircraft industry. A prime 
example is the Blended- Wing-Body (BWB) which 
has established a renaissance in the design of a family 
of all-new aircraft [4]. This revolutionary concept is 
enabling aerodynamic advances in all of the above 
areas, and then some. It presents challenges, yet 
offers significant opportunities, and as a result, a 5% 
reduction in DOC is within grasp. Suffice it to say 
that aerodynamics is not a sunset technology, but 
rather, it is as important today as it was a century 
ago; only the stakes have changed. 


2 Background 

The first and second authors embarked on a collab- 
orative effort which began in late 1994. Specifically, 
this study was to determine what was required to 
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obtain accurate drag results from the OVERFLOW 
code [5]. At that time, a typical simulation for a 
commercial transport configuration yielded a drag 
error of about 100 counts, where the total drag of 
the aircraft was nominally 300 counts. Results for 
the High-Speed Research (HSR) platforms were even 
worse; computed drag values were occasionally neg- 
ative. Clearly, that state-of-the-art was quite unac- 
ceptable. In fact, there were prominent members of 
the aerodynamics community at large who felt that 
accurate drag predictions from RANS-based CFD 
methods might never be accomplished. Nonetheless, 
a systematic study of gridding guidelines proved to 
be the key, and by mid 1995 the errors in computed 
absolute drag values for pre-buffet cruise conditions 
approached the 1-3% level. This collaboration, along 
with the NASA Advanced Subsonic Transport (AST) 
Program cooperative work with (then) McDonnell 
Douglas Aerospace- West [6], was the catalyst for the 
development of the 7-zone grid system for wing/body 
configurations with the signature collar grid at the 
wing/body juncture. These gridding guidelines were 
used as the basis for the DPW baseline grids, and will 
be discussed later in the paper. However, the end re- 
sult was that the size of a grid suitable for accurate 
drag prediction was nominally 4-5 times larger than 
that previously used. In addition, the number of iter- 
ations required for convergence on drag rather than 
pressures also jumped by a factor of 2-5, depending 
on the case. Hence, the cost of OVERFLOW sim- 
ulations for drag prediction increased by more than 
an order-of-magnitude relative to those used for the 
calculation of pressure distributions. 

In another collaborative effort which began in 
early 1995, the first two authors agreed that a paral- 
lel version of OVERFLOW for distributed process- 
ing was in order. The original parallel code was 
based on the Parallel Virtual Machine (PVM), and 
more recently, has been ported to use the Message- 
Passing Interface (MPI) [7]. For more than five 
years, the Aerodynamic Design group in Long Beach, 
CA has almost exclusively used a parallel version of 
OVERFLOW on distributed clusters for production 
overset-grid CFD simulations. A parallel-processing 
capability such as this was necessary for accurate 
CFD drag predictions to be economically feasible 
in an aircraft design environment [8]. Turn-around 
time was further improved with the addition of grid 
sequencing and full multigrid to accelerate solution 
convergence [9]. 

The combination of the above two collaborations 
had an immediate impact on the B7 17-200 design 
(previously called the MD-95). Here, several aerody- 
namic fairing designs for various juncture flows were 
evaluated using OVERFLOW. In all cases, the pre- 
dicted drag increments were later confirmed to be 


extremely accurate in wind-tunnel tests [10]. In one 
particular case, a pocket of separated flow was identi- 
fied in the numerical simulations just prior to a wind- 
tunnel entry. This prompted a flow-visualization 
run that confirmed the separation. Before the test 
was over, a fillet for this troublesome region was de- 
signed, fabricated with stereolithography, shipped to 
the wind tunnel, and tested. 

In 1996, the first author invited a team of NASA 
personnel from the Ames and Langley Research Cen- 
ters to participate in the aerodynamic design of the 
MD-XX trijet aircraft. The NASA group worked on- 
site and fully integrated within the MD-XX team. 
Due to the successes of CFD drag prediction on the 
MD-95 program, the MD-XX Design Office elevated 
the role and importance of CFD in the design envi- 
ronment. They did so by scheduling the freeze of the 
final loft lines of the cruise geometry several months 
prior to the first wind-tunnel test entry which would 
verify the aircraft’s aerodynamic performance. Al- 
though this program was later cancelled, it marked 
a dramatic change in the manner in which business is 
conducted in the design of an all-new aircraft. This 
philosophy lives on today in advanced programs such 
as the Blended- Wing-Body [4], [11]. 

The lessons learned in the above efforts have been 
augmented with subsequent studies on drag pre- 
diction conducted under various programs and ex- 
tended to other CFD methods such as CFL3D [12], 
SYN107 [13]-[14], and TLNS3D [15]. Within Boeing 
Phantom Works Long Beach, OVERFLOW remains 
the work-horse for complex transport configurations, 
CFL3D is heavily used on the BWB and re-entry ve- 
hicle programs, while SYN107 and TLNS3D round 
out the tool chest by providing aerodynamic shape 
optimization capabilities. 

Since 1995, while the errors in predicted absolute 
drag have stablized, the complexity of the configura- 
tions being analyzed has consistently increased. To- 
day, the size of an overset grid system for a complete 
B747-400 configuration (comprised of a cruise wing, 
fuselage, pylons, flow-through bifurcated fan and 
core cowls, winglets, vertical and horizontal tail com- 
ponents) is approximately 20 million nodes. Further- 
more, these simulations are being performed with the 
aircraft trimmed to specified center-of-gravity loca- 
tions. After correcting for excrescences, internal cowl 
drags, etc., comparisons with flight test data have 
confirmed that the numerically predicted absolute 
drag values are within the band of uncertainties. 

The next challenge for drag prediction is to im- 
prove the level of accuracy for post-buffet cruise con- 
ditions as well as for high-lift configurations. Im- 
proved high-lift drag prediction may become criti- 
cal to achieve the pending more-stringent environ- 
mental requirements on take-off and landing noise. 
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Unfortunately, current state-of-the-art RANS-based 
CFD methods cannot consistently predict accurate 
pressures for these flows. Until this is accomplished, 
there is little hope that accurate drag values will be 
attained here as well. An accomplishment of this 
magnitude will likely be possible only through the 
cumulative work of many collaborative efforts such 
as those aforementioned. 

The works noted above were conducted by a mul- 
titude of individuals, including the authors, under 
numerous collaborations between Boeing Phantom 
Works Long Beach and the NASA Ames and Lan- 
gley Research Centers. A subset of those who were 
involved are: Dan Bencze, Bob Biedron, Dick Camp- 
bell, William Chan, Roger Clark, Susan Cliff, Mark 
DeHaan, Lie-Mine Gea, Robb Gregg, James Hager, 
Ray Hicks, Rick Hooker, Dennis Jespersen, Steve 
Krist, Steve Mysko, Bob Narducci, Mike Olsen, Rick 
Potter, Stuart Rogers, Dino Roman, Tony Sclafani, 
Jeff Slotnick, Richard Wahls, and Mark Whitlock. 

3 DLR-F4 Geometry 

The case chosen for the DPW is the DLR-F4 
wing/body configuration [16]. Several factors were 
considered in the decision to use this geometry as 
the test-bed for the workshop. One factor was the 
availability of test data from multiple wind-tunnel 
facilities. Another was that this configuration is rep- 
resentative of current transonic transport aircraft. 

The general layout of the DLR-F4 is provided in 
Figure 1. This configuration is typical of a transonic 
aircraft designed to cruise at M = 0.75. The wing 
quarter-chord is swept 25° with a leading-edge sweep 
of 27.1° and an outboard trailing-edge sweep of 18.9°. 
The 9.5 aspect-ratio wing is rigged with a dihedral 
angle of 4.8°. Its planform is void of a leading-edge 
glove, yet includes a yehudi which extends to 40% 
semispan, completely unsweeping the trailing edge 
of the inboard wing. This planform is representa- 
tive of wings that accommodate retractable main 
landing gear. The airfoil sections are supercritical 
with thickness-to-chord ratios of 14.9% at the side- 
of-body, reducing to 12.2% outboard. The wing trail- 
ing edge has a blunt base of 0.5% local chord. Fig- 
ures 2-3 provide the thickness and camber distribu- 
tions, as well as the geometry, for the root and out- 
board airfoil sections, respectively. The wind-tunnel 
model has a wing semispan of 585.7mm, a mean 
aerodynamic chord of 141.2mm, a reference area of 
145,400mm 2 , and reference center at x = 504.9 mm. 
The fuselage length is 1,192mm. Its constant bar- 
rel section has a diameter of 148.42mm, which be- 
gins at x = 250mm and extends to x — 626mm. 
No special fillets are incorporated at the wing-body 


juncture, yielding a sharp corner everywhere on the 
intersection line. The DPW geometry also includes 
the aeroelastic twist of the wind-tunnel model under 
a loading which corresponds to the nominal cruise 
condition of M - 0.75 and Cl = 0.5, with a dy- 
namic pressure of q = 43, 434Pa. 

4 DPW Overset Grid 

The overset mesh generated for the DPW was based 
on the original process developed in 1995. The grid 
is comprised of 7 zones, 4 of which conform to the 
geometry and three box grids that transition the sys- 
tem to the farfield boundary. The 4 conforming grids 
define the volumes next to the fuselage, the wing- 
body juncture, the wing, and the wingtip. Two in- 
termediate boxes surround the fuselage and wing ge- 
ometries, while the remaining farfield box extends 
outward about 150 reference-chord lengths. 

For this exercise, the surface grids were con- 
structed using Gridgen-V13 [17] and are depicted 
in Figure 4. These surface grids were then ex- 
truded outward using HYPGEN [18] to generate the 
4 surface-abutting volume grids. Figure 5 shows a 
close-up of the wingtip grid. Figures 6-7 illustrate 
the collar grid at the wing-body juncture. 

The hole-cutting and fringe-point coupling steps 
were performed manually using GMAN [19]. The 
overlap and blanking of the meshes near the wing’s 
mid-chord are shown in Figure 8. 

On the wing surface, the chordwise spacing at both 
the leading and trailing edges is approximately 0.1% 
local chord. The trailing-edge base is defined with 5 
evenly-spaced points. The wake is represented with 
65 points in the streamwise direction. The span wise 
spacing is about 1% semispan at the root and 0.1% 
at the tip. On the fuselage nose and after-body, the 
maximum grid spacing is nominally 5mm. In the 
direction normal to the viscous surfaces, the first- 
layer spacing is about 0.001mm, which corresponds 
to — T Also in this direction, the maximum 
growth rate of the grid spacing is 1.24. Figure 9 
provides an itemization of the individual grid di- 
mensions, surface points, total grid points, and non- 
blanked real points. The complete grid system is 
comprised of 3,231,377 real points, with 54,445 re- 
siding on the viscous surfaces. 

The guidelines used to generate the DPW overset 
mesh purposely omitted two gridding rules; these will 
be discussed now for completeness. The first is re- 
lated to the manner in which OVERFLOW computes 
skin-friction drag. For this calculation to be second- 
order accurate, the first two layers of the grid nor- 
mal to viscous walls must be evenly spaced. While 
this rule was not strictly enforced, the spacing ra- 
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tios of the first two layers were fairly close to unity. 
The second rule is related to the grid resolution at a 
blunt trailing-edge base. Here, it is the first author’s 
standard practice to include a trailing-edge cap grid 
that wraps around the blunt base in a C-clamp fash- 
ion. This grid normally has half of the cells on the 
base and the remaining cells evenly split between the 
upper and lower surfaces. It normally extends about 
5-10% upstream of the trailing edge. However, in the 
case of the DLR-F4 wing where the gradients near 
the trailing edge are relatively benign, inclusion of a 
high-resolution cap grid is probably not required for 
accurate drag prediction. 

One final note. After the DPW was held, the 
Long Beach Aerodynamic Design group has finally 
transitioned from GMAN to PEGASUS-V5 [20] to 
provide semi-automatic hole-cutting and fringe-point 
coupling capabilities. There were several reasons 
for this lag, most were related to drag prediction. 
With that being said, the current parallel version of 
PEG ASUS- V5 has proved to be a very useful tool, 
and one now appropriate for drag prediction. Fur- 
ther, it will continue to improve through the on-going 
collaborations within the PEGASUS community. 


5 Wind-Tunnel Test Data 

The DLR-F4 model was tested in three European 
facilities: NLR-HST, ONERA-S2MA, and DRA-8x8. 
The repeatability of these facilities was on the order 
of ±5 counts. While the AGARD AR-303 report [16] 
presented this data, an unfortunate element of this 
documentation was that the drag coefficients were 
only tabulated to three decimal places. Hence, the 
archived public-domain data has an effective scatter 
band of ±10 counts. 

In an attempt to alleviate the uncertainty in- 
troduced by the truncation, the first author post- 
processed the public-domain data with two filters. 
The first data enhancement augmented the drag co- 
efficients with the coefficients of axial and normal 
forces. This process reduced the uncertainty of the 
tabulated values from the original 5 counts, down 
to as small as 0.8 counts, depending on the case. 
The second data enhancement was a careful digiti- 
zation of the drag polar figure in the AGARD re- 
port. This digitization was also checked to be consis- 
tent with the reduced uncertainty bands derived by 
the first enhancement process. The resulting wind- 
tunnel drag polars are provided in Figure 10. A more 
detailed explanation of these filters can be found in 
the Data Enhancement presentation on the DPW 
website [1]. 

The enhanced wind-tunnel test data for M = 0.75 
has been collapsed to a curve, by fitting a limited 


range of the data to an equation of the following 
form. 


Cd = Coo + 


Cl 


( 2 ) 


7r * e * AR 

Here, Coo and t are the free coefficients of the curve 
fit and the aspect ratio of the DLR-F4 wing is AR = 
9.437262. Including data from all three tests, but 
limited to the lifting range of 


0.2 < C L < 0.55, (3) 

a least-squares curve fit yields Coo = 0.018388 and 
e = 0.829146. Eqn (2) now becomes: 

Cl> = 0018388 + 24 ^ 55 ’ <“> 

and is applicable for the Cl range given in Eqn (3). 
The curve fit given by Eqn (4) and associated test 
data are provided in Figure 11. 

Since Case 1 of the DPW exercise is to compute 
the drag at M = 0.75, Ren = 3 AT and Cl - 0.5, it 
might be of interest to apply Eqn (4) to this condi- 
tion, which yields; 


C,, = 0^8388+ 3^ = 0.0S8M. < 5 > 

The drag polars presented in the next section use 
the enhanced wind-tunnel data, rather than data 
taken directly from the AGARD AR-303 report. 


6 Results 

Included in this section are comprehensive sets of 
OVERFLOW solutions which were performed on the 
DPW baseline overset mesh. Also included is a lim- 
ited set of data generated by CFL3D on the same 
overset mesh for additional comparison. A more 
complete set of CFL3D solutions for the DPW base- 
line 1-to-l multiblock mesh is documented by the 
third author in Reference [21]. 

For the DPW exercise, the authors purposely did 
not coordinate with each other in an attempt to 
obtain independent results on the DLR-F4 config- 
uration. In spite of this, the first two authors ran 
OVERFLOW with essentially the same set of critical 
input parameters. Both used central difference scalar 
dissipation with the Spalart-Allmaras (SA) turbu- 
lence model. However, slightly different versions of 
OVERFLOW were run, different computer platforms 
were utilized at different levels of precision, and dif- 
ferent methods to converge on lift were employed. It 
is good to report that these differences yielded no no- 
ticeable variations in the computed forces, moments 
or pressures. 

The first author ran full convergence on all so- 
lutions, starting each solution from a uniform flow 
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at freestream conditions. Full multigrid acceleration 
was used with 150 iterations in both the coarse and 
intermediate meshes, and 3, 000 iterations in the fine 
mesh. All solutions were run by specifying an angle- 
of-attack, even when a specific lifting condition was 
desired. All computations were performed using 64- 
bit precision. Version 1.8m was run in parallel using 
MPI on a cluster of 6 Hewlett-Packard C3610 work- 
stations, each with 2 GB of RAM, and connected 
with a switched 100BaseT ethernet. Each solution 
required about 13 hours of wall-clock time. Alpha 
sweeps were conducted at 10 Mach numbers, then 
these data were interpolated on Cl to derive the re- 
quired alphas for C L s of 0.3, 0.4, 0.5, and 0.6. In all, 
a total of 53 OVERFLOW solutions were performed 
to obtain drag polars at the 10 Mach numbers. To 
define the drag-rise curves, the drag polars were in- 
terpolated on C\ to obtain the corresponding Cl 
values. 

The second author and colleagues ran some cases 
with a fixed angle-of-attack and some with a specified 
lift-coefficient. All computations were performed us- 
ing 32-bit precision. Version 1.8s was run in parallel 
on three different computer platforms, an SGI Oc- 
tane with 2 processors, an SGI Origin using 8 proces- 
sors, and a cluster of 6 Compaq XP-1000 machines. 
Wall-clock timings for these systems were 18.5 hours, 
7.5 hours, and 6 hours per 1, 000 fine-mesh iterations, 
respectively. 

For a limited set of conditions on the overset mesh, 
the third author ran CFL3D with 3 r ^-order upwind 
differencing with Roe flux difference splitting and us- 
ing the SA model. These solutions were converged 
4,800 multigrid iterations on the fine grid only, and 
were run by specifying an angle-of-attack. A non- 
dedicated SGI Origin was used with one processor, 
requiring approximately 270 hours of wall-clock time 
per solution. 

Figures 12-13 illustrate the OVERFLOW conver- 
gence histories of lift and drag, respectively, for a 
freestream condition of M = 0.75, Ren = 3 M, and 
a = -1°, which yields Cl = 0.409 and Co = 263.3 
counts. These forces have essentially converged by 
2,000 iterations on the fine mesh. Note the scale 
on these figures; lift increments are 0.002 and drag 
increments are a count. 

Figures 14-15 provide a similar set of OVER- 
FLOW convergence histories, however, these corre- 
spond to the cruise lifting condition of C L - 0.5 
with corresponding C D - 295.6 counts. Unlike the 
previous condition, these forces continue to oscillate 
through all 3, 000 iterations, albeit at very small am- 
plitudes. The cause of these fluctuations is a small 
pocket of separated flow that appears near the trail- 
ing edge of the wing-body intersection; see Figure 16. 
Figure 17 provides a comparison of pressure distri- 


butions between OVERFLOW and test data at the 
cruise design point of M = 0.75 and C L - 0.5. Note 
that the leading-edge peaks are missed because of 
the difference in alphas between the numerical sim- 
ulation and the tests. Also, the isobars in this figure 
hint to the small pocket of flow separation on the 
upper-surface near the root trailing edge. 

An important aspect of comparing results from 
OVERFLOW and CFL3D on the overset mesh is to 
estimate variation due to choice of CFD code. Fig- 
ure 18 shows a comparison of pressure distributions 
at M = 0.75 and a = 0°. At this a-matched con- 
dition, results are very close and computed leading- 
edge peaks better match the test data. The lift, drag 
and moment comparisons are provided in Table 1. 
Additional documentation on variations due to CFD 
method can be found in Reference [22]. 



C L 

Cl ( counts ) 

Cm 

OVERFLOW 

0.532 

309.6 

-0.1614 

CFL3D 

0.535 

313.0 

-0.1653 


Table 1: OVERFLOW and CFL3D Comparison at 
M ~ 0.75, Ren ~ 3 M, and a = 0°. 


Case 1 of the DPW asked that all participants 
compute a solution for the cruise design point of 
M — 0.75, Ren — 3 Af, and Cl — 0.5. The sec- 
ond author and two colleagues independently ran 
this case and included solutions with 3rd-order R.oe 
upwind as well. Combined with the first author’s 
data, these results are provided in Table 2. 


Alpha (deg) 

Cl ( counts ) 

Type 

-0.260 

295.6 

Central 

-0.258 

295.8 

Central 

-0.257 

295.9 

Central 

-0.254 

296.0 

Central 

-0.241 

292.6 

Roe 

-0.238 

292.9 

Roe 


Table 2: DPW Case 1 OVERFLOW Results. 

This table shows that the independent assessments 
of drag using a consistent differencing scheme were 
within ±0.2 counts of the average. Correcting for 
the variation of Cl (which occurs in the 4 th decimal 
place) collapses the variations in drag to an insignif- 
icant number. The variation between central differ- 
encing and Roe upwinding is ~ 3 counts, with Roe 
yielding the lower drag. 

While a CFL3D solution for Cl = 0.5 was not 
computed on the overset mesh, Table 3 contains data 
which stradle this condition, and include interpo- 
lated values for a (based on Cl) and Cl (based on 
Cl). This table shows that the interpolated CFL3D 
result yields a drag level approximately 2-5 counts 
higher than OVERFLOW at the design point. 


Vassberg, Buning & Rumsey, AIAA Paper 2002-0840, Reno, NV 


6 of 23 





Alpha (deg) 

Co ( counts ) 

Cl 

- 1.000 

266.0 

0.416 

-0.294 

298.0 

0.500 

0.000 

313.0 

0.535 


Table 3: CFL3D Results. 


Figure 19 provides the OVERFLOW computed 
drag polar for M = 0.75 and includes the wind- 
tunnel data for reference. The data in this figure are 
inconsistent in that the numerical simulations were 
performed assuming fully-turbulent flow, while the 
tests allowed laminar runs of 5-15% chord on the 
wing’s upper surface and 25% chord on the wing’s 
lower surface. Figure 20 shows the transitions pat- 
tern used in the tests. 

To assess the impact on drag caused by this dif- 
ference, FL022 [23]- [24] was used to generate drag 
polars of fully-turbulent flows and flows tripped with 
the pattern of Figure 20. These results are shown 
in Figure 21. Also shown in this figure are least- 
squares curve fits of the polars. When differenced, 
these curve fits yield a shift in drag of 

Cd shift — 13.7 — 2.3 * Ci counts. (6) 
Hence at Cl = 0.5, 

C Ds hijt{ 0.5) = 13.7- 2.3 *0.5 2 

= 13.1 counts . (7) 

Applying the correction of Eqn (6) to the fully- 
turbulent OVERFLOW and CFL3D results yields 
the drag polars depicted in Figure 22. Once the effect 
of transition is taken into consideration, the numer- 
ically predicted results lie within the scatter of the 
experimental data. 

If the cruise point is reviewed, the drag predictions 
with the correction of Eqn (7) become: 

Cd overflow corrected — - 0.02825, 

C D cfl3d—v6 corrected 0.02849. (8) 

Comparing Eqns (8) with Eqn (5) shows only a 1% 
difference in drag levels between CFD and experi- 
ment. 

During the DPW, Schwamborn and Sutcliffe pre- 
sented a result using the DLR-TAU code which also 
addressed this CFD-to-Test difference. Their presen- 
tation is available on the DPW T website [1]. On page 
5, the data labeled Case S and Transition indicate 
that at Cl = 0.5, the correction is: 

Cd j/u/t(0.5) ~ 12 counts. (9) 

Eqns (7) & (9) are very consistent with each other 

and are independent estimates of the effects on drag 
of the laminar runs vs. fully-turbulent flows. 


Figure 23 illustrates the drag-rise curves as pre- 
dicted by the OVERFLOW solutions. From bottom- 
to-top, the four curves of this figure correspond to 
Cl = 0.3, 0.4, 0.5, and 0.6, respectively. Also in- 
cluded in this figure is a cross-plot curve, depicted by 
the dotted line, which illustrates the drag-divergence 
Mach number for the DLR-F4. The definition of 
Mdd is taken to occur when the drag-rise slope is 
dCo/dM = 0.05. Using these results, the drag- 
divergence boundary is determined and is shown in 
Figure 24. 

The above definition of Mdd is motivated par- 
tially by the Breguet- Range Eqn (1), and by eco- 
nomic forces (related to block times) that push the 
operating Mach number upward towards the 99% 
Long-Range-Cruise (LRC) point. From Figure 24, it 
appears that the DLR-F4 wing is capable of cruis- 
ing in excess of M ~ 0.775 at the lifting condi- 
tion of Cl — 0.5. However, these estimates have 
all been made at a wind-tunnel Reynolds number, 
rather than at flight. 

7 Conclusions 

The importance of drag prediction within an aircraft 
design environment is reviewed. A seven-year ret- 
rospective of collaborations is given which illustrates 
some of the steps taken by the authors and colleagues 
in pursuit of accurate drag predictions. The suc- 
cesses of this body of work have had a significant 
impact on the manner in which all-new aircraft de- 
signs are approached. 

The DLR-F4 wing/body configuration has been 
analyzed by OVERFLOW and CFL3D using an over- 
set mesh. This study focused on the prediction of 
drag as a function of lift and Mach number for a 
wind-tunnel Reynolds number of 3 M based on ref- 
erence chord. For pre-buffet conditions, results pre- 
sented herein show that the numerical drag predic- 
tions are within 1% of the averaged wind-tunnel 
data, and the OVERFLOW and CFL3D drag pre- 
dictions are within 1-2% of each other. This level of 
uncertainty is comparable to that of the wind-tunnel 
data itself. 

Comprehensive sets of OVERFLOW simulations 
were independently performed by several users, on 
a variety of computational platforms. These solu- 
tions spanned 10 Mach numbers from 0,5 to 0.82. 
By coincidence, most of these solutions used cen- 
tral differencing, scalar dissipation and the Spalart- 
Allmaras turbulence model. However, a variety of 
parallel computational platforms were used, a mix- 
ture of single and double precision simulations were 
performed, and slightly different versions of OVER- 
FLOW were applied. It is good to say that the vari- 
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ation of these results are negligible. In addition, a 
few computations were performed with Roe upwind- 
ing. The difference in drag levels between the two 
differencing stencils is on the order of 1%. 

The Drag Prediction Workshop called for the CFD 
calculations to be performed fully turbulent. The 
available wind-tunnel data, however, was tripped at 
5-15% on the wing upper surface and 25% on the 
lower surface. Corrections for the CFD-to-Test dif- 
ferences were estimated using FL022, running semi- 
complete polars for both flows. The study derived a 
correction of about 13 counts at the design point. 
This is comparable to a similar and independent 
study based on the DLR-TAU code. Drag polars 
based on OVERFLOW were constructed with the 
CFD-to-Test corrections and compared with wind- 
tunnel test data. These comparisons show that the 
predicted polars fall within the scatter band of the 
test data, at least for pre-buffet conditions. 

Four drag-rise curves were constructed from the 
OVERFLOW solutions. Using a slope definition for 
drag-divergence Mach number, the Mdd boundary 
for the DLR-F4 wing/body was constructed. 

In future workshops, a grid-resolution study 
should be included by providing a series of 
parametrically-consistent meshes of varying resolu- 
tion. Post-processing the results from such a se- 
quence using Richardson extrapolation will provide 
further insight into the resolution required for a de- 
sired level of accuracy, at least for the case being 
investigated. 

Acknowledgment 

The authors acknowledge The Boeing Company and 
The National Aeronautics and Space Administration 
for their support in the Drag Prediction Workshop, 
and for providing a working environment that not 
only allows, but encourages the numerous collabo- 
rative studies that have and will continue to occur 
between the two institutions. 

The second and third authors would like to ac- 
knowledge S. Melissa Rivers, Joseph H. Morrison, 
and Robert T. Biedron of the NASA Langley Re- 
search Center, for their contributions to the authors’ 
respective Drag Prediction Workshop presentations. 

References 

[1] AIAA CFD Drag Prediction Workshop Website. 
http://www.aiaa.org/tc/apa/dragpredworkshop 
/dpw.html, June 2001. 

[2] D. W. Levy, J. C. Vassberg, R. A. Wahls, 
T. Zickuhr, S. Agrawal, S. Pirzadeh, and M. J. 


Hemsch. Summary of data from the first AIAA 
CFD Drag Prediction Workshop. AIAA pa- 
per 2002-0841 , Reno, NV, January 2002. 

[3] J. C. Vassberg and R. D. Gregg. Overview of 
aerodynamic design for transport aircraft. Pre- 
sentation , First M.I.T. Conference on Compu- 
tational Fluid and Solid Mechanics, Cambridge, 
MA, June 2001. 

[4] R. H. Liebeck. Design of the Blended-Wing- 
Body subsonic transport. Wright Brothers Lec- 
ture, AIAA paper 2002-0002 , Reno, NV, Jan- 
uary 2002. 

[5] P. G. Buning, D. C. Jespersen, T. H. Pulliam, 
G. H. Klopfer, W. M. Chan, J. P. Slotnick, S. E. 
Krist, and K. J. Renze. OVERFLOW User’s 
Manual, Version 1.8L. Technical report, NASA, 
July 1999. 

[6] L. M. Gea, N. D. Halsey, G. A. Intemann, 
and P. G. Buning. Applications of the 3-D 
Navier-Stokes code OVERFLOW for analyz- 
ing propulsion-airframe integration related is- 
sues on subsonic transports. In Proceedings of 
19th Congress of the International Council of 
the Aeronautical Sciences (ICAS 94), number 
ICAS-94-3.7.4, pages 2420-2435, Anaheim, CA, 
September 1994. 

[7] D. C. Jespersen. Parallelism and OVER- 
FLOW. NAS Technical Report NAS-98-013, 
NASA Ames Research Center, October 1998. 
http://www.nas.nasa.gov/Research/Reports 
/Techreports/1998/PDF/nas-98-013.pdf. 

[8] J. Conlon. OVERFLOW code empowers Com- 
putational Fluid Dynamics. InSights Issue 5, 
NASA High-Performance Computing and Com- 
munications Program, Moffett Field, CA, April 
1998. 

[9] D. C. Jespersen, T. H. Pulliam, and P. G. Bun- 
ing. Recent enhancements to OVERFLOW. 
AIAA paper 97-0644 , AIAA 35th Aerospace Sci- 
ences Meeting, Reno, NV, January 1997. 

[10] J. C. Vassberg. CFD-Based Design. Presenta- 
tion, Stanford University, Palo Alto, CA, Febru- 
ary 2000. Graduate Seminar Series. 

[11] D. L. Roman, J. B. Allen, and R. H. Liebeck. 
Aerodynamic design challenges of the Blended- 
Wing-Body subsonic transport. AIAA pa- 
per 2000-4335 , Denver, CO, August 2000. 

[12] S. L. Krist, R. T. Biedron, and C. L. Rum- 
sey. CFL3D User’s Manual. NASA-TM 1998- 
208444, NASA, June 1998. Version 5.0. 


Vassberg , Buning & Rumsey, AIAA Paper 2002-0840 , Reno, NV 


8 of 23 



[13] A. Jameson and J. C. Vassberg. Computational 
fluid dynamics for aerodynamic design: Its cur- 
rent and future impact. AIAA paper 2001-0538 , 
Reno, NV, January 2001. 

[14] J. C. Vassberg and A. Jameson. Aerodynamic 
shape optimization of a Reno race plane. Inti 
Journal of Vehicle Design , 28(4) :318— 338, 2002. 
Special Issue on: Design Sensitivity and Opti- 
mization. 

[15] G. Kuruvilla, R. P. Narducci, and S. Agrawal. 
Development and application of TLNS3D- 
Adjoint: A practical tool for aerodynamic shape 
optimization. AIAA paper 2001-2400, Anaheim, 
CA, June 2001. 

[16] G. Redeker. DLR-F4 wing-body configuration. 
In A Selection of Experimental Test Cases for 
the Validation of CFD Codes , number Alt-303, 
pages B4.l-B4.21. AGARD, August 1994. 

[17] Gridgen user manual version 13. Technical re- 
port, Pointwise, 1998. 

[18] W. M. Chan, I. T. Chiu, and P. G. Buning. 
User’s manual for the HYPGEN hyperbolic grid 
generator and the HGUI graphical user inter- 
face. NASA-TM 108791, NASA, October 1993. 

[19] T. D. Gatzke, W. F. LaBozzetta, G. P. Frin- 

frock, J. A. Johnson, and W. W. Romer. 
MACGS: A zonal grid generation system for 

complex aero-propulsion configurations. AIAA 
Paper 91-2156, June 1991. 

[20] N. E. Shus, W. E. Dietz, S. M. Nash, M. D. 
Baker, and S. E. Rogers. PEGASUS user’s 
manual version 5.1c. Technical report, MICRO 
CRAFT, July 2000. 

[21] C. L. Rumsey and R. T. Biedron. Compu- 
tation of flow over a drag prediction work- 
shop wing/body transport configuration using 
CFL3D. NASA-TM 2001-211262, NASA, De- 
cember 2001. 


[24] P. A. Henne and R. M. Hicks. Wing analysis 
using a transonic potential flow computational 
method. NASA-TM 78464 , July 1978. 


[22] C. L. Rumsey, D. O. Allison, R. T. Biedron, 
P. G. Buning, T. G. Gainer, J. H. Morrison, 
S. M. Rivers, S. J. Mysko, and D. P. Witkowski. 
CFD sensitivity analysis of a modern civil trans- 
port near buffet-onset conditions. NASA-TM 
2001-211263, NASA, December 2001. 

[23] A. Jameson. Transonic potential flow calcula- 
tions using conservative form. In Proceedings of 
AIAA 2nd Computational Fluid Dynamics Con- 
ference, pages 148-161, June 1975. 


Vassberg , Buning & Rumsey, AIAA Paper 2002-0840, Reno , NV 


9 of 23 



SZ'ltll 



Figure 1: General Layout of the DLR-F4 Model. 
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Figure 2: Root Airfoil Section of the DLR-F4 Model. 
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Figure 6: Collar Grid Volume 



Figure 7: Collar Grid on Fuselage Near Wing Trailing Edge 
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Figure 8: Field Grids Near Wing Mid-Chord Location 


Total viscous-surface points: 54445 
Total grid points: 3727462 
Total non-blanked grid points: 3231377 
Grid points across the TE base 5 


Farfield boundary is a box -150 chord lengths away. 


grid 

id 

jd 

kd 

surfpts 

gridpts 

%pts 

realpts 

%real 

description 

1 

49 

273 

49 

13377 

655473 

17.6 

562499 

17.4 

fuselage 

2 

385 

65 

49 

22209 

1226225 

32.9 

1092327 

33.8 

collar 

3 

385 

62 

49 

15934 

1169630 

31.4 

1038606 

32.1 

wing 

4 

25 

141 

49 

2925 

172725 

4.6 

135219 

4.2 

wingtip 

5 

121 

22 

43 

0 

114466 

3.1 

76840 

2 . 4 

fuse_box 

6 

67 

74 

46 

0 

228068 

6.1 

168408 

5.2 

wing_box 

7 

75 

39 

55 

0 

160875 

4.3 

157468 

4.9 

global_box 


Totals: 54445 3727462 3231377 


Figure 9: Statistics of the DPW Baseline Overset Grid. 
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DLR-F4 WING / BODY GEOMETRY 
AGARD Report 303 - Figure 12 

Mach = 0.75 , Re = 3 Million 



COEFFICIENT OF DRAG 
Digitized Data 

re 10: Enhanced Test Data M = 0.75 Drag Polar for the DLR-F4. 


DLR-F4 WING / BODY GEOMETRY 
AGARD Report 303 Data 


Mach = 0.75 , Re = 3 Million 



COEFFICIENT OF DRAG 
Digitized Data 


Figure 11: Least-Squares Fit of the DLR-F4 Enhanced Test Data. 
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DLR-F4 WING/BODY CONFIGURATION 
Baseline Over-Set Grid 

M = 0.75 , REN = 3M , Alpha* *1.00 , CL = 0.40864 , CD = 0.026332 



FINE-MESH ITERATIONS 


Figure 12: Convergence History of Lift at CL— 0.409. 

DLR-F4 WING/BODY CONFIGURATION 
Baseline Over-Set Grid 


M = 0.75 , REN = 3M , Alpha = -1.00 , CL = 0.40864 , CD = 0.026332 



FINE-MESH ITERATIONS 

Figure 13: Convergence History of Drag at CL=0.409. 
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DLR-F4 WING/BODY CONFIGURATION 
Baseline Over-Set Grid 



DLR-F4 WING/BODY CONFIGURATION 
Baseline Over-Set Grid 



Figure 15: Convergence History of Drag at CL— 0.500. 
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Figure 16: Juncture-Flow Separation at Design Point. 
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DLR-F4 WING/BODY CONFIGURATION 
CFD-TO-TEST CORRECTIONS 


M = 0.75 , REN = 3M 
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Figure 21: FL022 Drag Polars with and without Transition Trips. 

DLR-F4 WING/BODY CONFIGURATION 
Baseline Over-Set Grid 


M = 0.75 . REN = 3M 



DRAG COEFFICIENT 

Figure 22: OVERFLOW Drag Polars with and without Transition Corrections. 
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